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The Very Large Area gamma-ray Space Telescope (VLAST) is a mission concept proposed to detect gamma- 
ray photons through both Compton scattering and electron-positron pair production mechanisms, thus enabling 
the detection of photons with energies ranging from MeV to TeV. This project aims to conduct a comprehensive 
survey of the gamma-ray sky from a low-Earth orbit using an anti-coincidence detector, a tracker detector 
that also serves as a low-energy calorimeter, and a high-energy imaging calorimeter. We developed a Monte 
Carlo simulation application of the detector using the GEANT4 toolkit to evaluate the instrument performance, 
including the effective area, angular resolution, and energy resolution, and explored specific optimizations of 
the detector configuration. Our simulation-based analysis indicates that the current design of the VLAST is 
physically feasible, with an acceptance above 10 m? sr which is four times larger than that of the Fermi-LAT, 
an energy resolution better than 2% at 10 GeV, and an angular resolution better than 0.2° at 10 GeV. The 
VLAST project promises to make significant contributions to the field of gamma ray astronomy and enhance 


our understanding of the cosmos. 
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I. INTRODUCTION 


Gamma-ray astrophysics is an exciting field in astronom- 
ical science that has received significant attention. Detect- 
ing cosmic gamma-ray emissions in the energy range from 
MeV to GeV is difficult using ground-based telescopes and 
is preferably performed in space. Gamma ray telescopes 
covering this energy range can be roughly divided into two 
categories based on the detection principle: pair-production 
telescopes and Compton-scattering telescopes. For pair- 
production telescopes, the OSO-3 [1] provided the first confir- 
mation that the detection of gamma-rays is feasible for com- 
plex backgrounds of charged particles. Breakthrough discov- 
eries of high-energy gamma ray observations were carried 
out using SAS-2 [2] and COS-B [3] missions in the 1970s. 
In the 1990s, EGRET made significant progress in surveying 
the gamma-ray sky above 50MeV, leading to the discovery 
of numerous high-energy gamma-ray sources [4]. The Fermi 
LAT gamma-ray space telescope [5], launched in 2008, has 
been highly successful in this field for over a decade, iden- 
tifying more than 6000 gamma-ray sources in its fourth cat- 
alog [6]. However, owing to the limited acceptance and an- 
gular resolution of detectors, nearly one-third of the sources 
remain unidentified. The GAMMA-400 space mission, in- 
stalled on the Russian Navigator space platform, is currently 
under preparation. With its excellent energy resolution and 
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unprecedented angular resolution above 30 GeV compared to 
other space-based systems, it has the potential to unlock new 
insights in this field. However, its effective area is limited 
to 4000 cm?, owing to the size of the detector [7]. The pi- 
oneering Compton-scattering telescope for opening the MeV 
gamma-ray astronomical window was COMPTEL [8]. An- 
other small detector onboard the Chinese space station, PO- 
LAR [9] is dedicated to measuring the polarization of MeV 
gamma rays through Compton scattering. The COSI [10] 
project, which is funded by NASA’s Small Explorer program, 
is scheduled for launch in 2026 and features an exceptional 
energy resolution. The effective area of COSI remains small 
and perhaps does not meet the requirements of a powerful 
detector for MeV time-domain astronomy. Recent advance- 
ments in detection technology (semiconductors, scintillators, 
and time-projection chambers) have sparked growing interest 
in the MeV energy band. Consequently, several space-based 
gamma-ray missions have been proposed in recent years, such 
as PANGU [11], AMEGO [12], ee ASTROGAM [13], AdEPT 
[14], and GECCO [15], MAST [16], GRAMS [17], XGIS- 
THESEUS [18], Crystal Eye [19], and MASS [20]. With 
ongoing developments in detection technology and increas- 
ing scientific demand, there is a pressing need for gamma-ray 
telescopes with enhanced sensitivity. 


We propose the Very Large Area gamma-ray Space Tele- 
scope (VLAST) [21-24] with a significantly larger effective 
area, four times larger than Fermi-LAT. VLAST primarily 
consists of an Anti-Coincidence Detector (ACD), a Silicon 
Tracker, a low energy gamma-ray detector (STED), and a 
high-energy imaging calorimeter (HEIC), combining the ad- 
vantages of DAMPE [25] and APT [26]]. The STED design 
can measure both low-energy Compton scattering and high- 
energy pair production events, which change the traditional 
tungsten plate to a thin cesium iodide (CsI) tile, allowing 
model-independent control of the systematic detector uncer- 
tainties. VLAST can survey the gamma-ray sky from a low 
Earth orbit in the energy band from approximately 0.1 MeV 


to more than 1 TeV. 

The key scientific goals of the VLAST include (i) searching 
for dark matter signatures in galaxy cores, galaxy disks, and 
dwarf galaxies [27-31], (ii) monitoring special gamma-ray 
sources over time, such as Active Galactic Nuclei, gamma-ray 
bursts, millisecond pulsars, and supernovas [32-37], (iii) un- 
derstanding the origin and transportation of cosmic rays [38- 
40], (iv) using extragalactic diffuse gamma-rays and gamma- 
ray horizons to study the evolution of the universe [41-44], 
(v) testing fundamental physical laws, such as the Lorentz in- 
variance and equivalence principle [45-48], and so on. 

In this study, we evaluated the performance and optimized 
the design of VLAST using simulations. The simulation 
framework was developed based on the GEANT4 [49]. For 
event reconstruction, different energy or trajectory recon- 
struction algorithms were used to analyze the Compton scat- 
tering and pair production events. The key performance pa- 
rameters of VLAST were obtained based on detailed MC sim- 
ulations, such as acceptance, effective area, angular resolu- 
tion, energy resolution, and e-p discrimination. Furthermore, 
some specific design parameters were optimized, such as the 
threshold of ACD, the size of ACD, the width of the silicon 
strip pitch, and the number of CsI layers in the STED. 


II. BASELINE DESIGN OF VLAST 


VLAST follows the foundational principles and structures 
of earlier gamma detectors. The detector included three pri- 
mary subdetectors: ACD, STED, and HEIC. The ACD plays 
a critical role in rejecting the charged particle background 
and reducing backsplashes from high-energy events. STED, 
which shares design similarities with APT, replaces the tung- 
sten plate in the tracker detector with thin CsI tiles. This alter- 
ation allows for the measurement of Compton events and pair 
production events, thus enabling the tracking of their trajec- 
tories to reconstruct the direction of the incident gamma rays. 
Low-energy MeV photons are primarily measured through 
Compton scattering with electrons to determine their energy 
and direction of incidence. The HEIC measures the energy 
of incident particles and images the profile of the electromag- 
netic or hadron shower of high-energy particles, which is used 
to discriminate between electrons and protons and provides 
an estimated direction of the incident particle. Fig. 1 presents 
a schematic of the VLAST detector. 


A. Anti-Coincidence Detector 


The primary function of the ACD is to minimize the back- 
splash effect caused by high-energy photons and to reject the 
background of the charged particles. As the first barrier for 
rejecting charged particles, it is closely wrapped around the 
STED. The ACD is composed of 448 plastic scintillator de- 
tector tiles (20cm x20 cm x1 cm). Plastic scintillators have a 
high detection efficiency for charged particles and a low de- 
tection efficiency for photons. Therefore, they are often used 
as the ACD for gamma ray detection. To prevent particles 


from passing through the gap between the ACD without being 
identified, the gaps between the tiles are covered by flexible 
scintillating fiber ribbons. 

The backsplash effect was first discovered by the EGRET 
team. When a high-energy photon enters the calorimeter, it 
generates numerous secondary low-energy photons (ranging 
from 0.1 MeV to a few MeV), some of which may travel in 
the opposing direction through the STED and reach the ACD. 
This scattering causes the electrons to recoil,thus generating 
an electrical signal in the ACD. The intensity of this effect 
depends on many factors, such as the incident particle en- 
ergy, calorimeter thickness, and material composition. The 
intensity of the backsplash effect increases with the energy of 
the incident particles and the radiation length of the calorime- 
ter. To mitigate this effect, the Fermi-LAT detector divides 
the ACD into smaller units by utilizing the discrepancy be- 
tween the piece in the incident direction and that in the recoil 
to discriminate the recoil event. This modular design effec- 
tively suppresses the misinterpretation of photons attributed 
to the backsplash effect of high-energy photons. Following 
this scheme, the ACD in the VLAST design enables the dis- 
missal of fired tiles that are not aligned with the direction 
of photon incidence, thereby significantly mitigating the in- 
fluence of the backsplash effect. The backsplash effect in 
VLAST is shown in Fig. 4. 

The photon flux is lower than that of protons by a factor of 
10° and lower than that of electrons by approximately 10° in 
the observational energy band of VLAST. To detect photons 
in the complex cosmic-ray background of charged particles, 
VLAST must have an outstanding ability to distinguish be- 
tween photons and charged particles. When photons trans- 
form into electron-positron pairs in the STED, they can be 
discriminated from protons based on the shower pattern pro- 
duced by the electrons-positrons and protons in the HEIC. As 
described in Section IV E, HEIC has good electron-to-proton 
(e-p) discrimination capability, which significantly reduces 
the proton background. However, the discrimination between 
photons and electrons depends only on the ACD, which re- 
quires a minimum rejection fraction of 0.999 for charged par- 
ticles. Therefore, ACD plays a critical role in identifying pho- 
tons and distinguishing them from charged particles in the 
VLAST detector. 


B. Silicon Tracker and low Energy gamma-ray Detector 


The STED of VLAST is designed to fulfill three primary 
functions. First, it reconstructs particle trajectories with an 
accuracy exceeding 120 um for most incident particles. Sec- 
ond, it determines the charge of cosmic rays. Lastly, it con- 
verts the incoming photons into electron-positron pairs and 
detects both the photons and electrons resulting from Comp- 
ton scattering. It consists of eight layers of silicon trackers 
and low-energy gamma ray detector modules. Each module 
consists of a CsI layer on the top and a silicon strip module 
at the bottom. The silicon strip module contains two layers 
of silicon strip detectors. The single-sided silicon strips are 
arranged in the x- and y-directions with a pitch size of 120 
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Fig. 1. Schematic of the VLAST detector. The ACD covers the whole detector except the bottom. The STED has a total of 8 layers, each 
containing one sub-layer of CsI crystals connected to wavelength-shifting fibers and two sub-layers of silicon strip modules. The HEIC is 
placed below the STED, which contains four layers of orthogonally arranged Bismuth Germanium Oxide (BGO) crystals. The bottom right 
panel shows the two classes of events detectable by VLAST, the Compton scattering events and the electron-positron pair production events. 


um, defined as the distance between the centers of the ad- 
jacent strips. The top six layers of the 8-layer CsI have a 
thickness of 2 mm each, and the bottom two layers have a 
thickness of 4 mm each, resulting in a total radiation length 
of approximately 1 Xo. Each layer of CsI is assembled by 
joining numerous CsI (Na) scintillating crystal square panels, 
with each panel side measuring approximately 200 mm. The 
upper and lower surfaces of each layer of CsI are tightly cou- 
pled to two layers of square wavelength-shifting (WLS) fibers 
arranged along the x- and y-axes, each with a cross-sectional 
side length of 2 mm. The WLS fibers absorb the blue scintil- 
lation light emitted from the CsI and transmit a fraction of the 
re-emitted red light to the SiPMs connected at their ends. The 
readout signals from the SiPMs can be used to infer the en- 
ergy deposition and x — y coordinates of particle interactions 
in the CsI panels. 


There are two main types of interactions between gamma- 
rays in the STED: Compton scattering and pair produc- 


tion. The dominant interaction depends on the energy of the 
gamma-rays. Compton scattering is predominant at energies 
below a few tens of MeV, whereas electron-positron pair pro- 
duction dominates at higher energies. Different types of de- 
tectors have been designed to capture gamma-rays in vari- 
ous energy bands. For example, COMPTEL [8], Fermi-LAT 
[5], and AGILE [50] focused on either the low-or the high- 
energy range. In contrast, VLAST aims to simultaneously 
detect gamma rays in the energy band simultaneously from 
0.1 MeV to 1 TeV, encompassing both interactions. Conse- 
quently, the gamma rays undergoing these two distinct in- 
teractions require different detection methods. The detection 
principle of VLAST is illustrated in the lower-right panel of 
Fig. 1. 


For pair production events, the conceptual design of 
VLAST is similar to that of Fermi-LAT and AGILE, but is 
optimized for the lower energy band. CsI serves a dual role as 
both a positron-electron pair converter and detector for mea- 


suring the position and low-energy deposition of interactions. 

This design ensures a high conversion efficiency of gamma 
rays into electron-positron pairs. The CsI low-energy 
calorimeter does not reduce the multiple scattering effects of 
pair events. However, it can provide the deposited energy and 
position of the fired tiles better than tungsten plates, allowing 
for a more accurate assessment of energy and multiple scat- 
tering. Consequently, this design improves the energy and 
angular resolution at low energies compared to Fermi-LAT 
and AGILE. 

Detection of a Compton scattering event is challenging be- 
cause the scattered photons carry away a significant portion of 
the momentum of the primary particle, especially when com- 
pared to pair production processes. Therefore, VLAST must 
be capable of simultaneously measuring two photons. Upon 
the entry of an incident photon into the detector, Compton 
scattering occurs within one layer of the tracker, giving rise 
to secondary photons and electrons. The photon is scattered 
with an electron in the detector and transfers a fraction of its 
energy (E1) to an electron. The scattered photon retains its 
remaining energy (£2) and may interact with the CsI of the 
lower-energy detector behind it or enter the HEIC. The angle 
0 between the photon and the electron after Compton scatter- 
ing is calculated as follows: 


1 1 
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where mec? is the mass-equivalent energy of electron. 

If only the angle of Compton scattering can be recon- 
structed, the incident direction of the photon can only be lo- 
cated on a circle in the sky, which is referred to as an “event 
circle”. These events are known as “untracked” events. The 
width of the circle is related to the accuracy of the detector in 
measuring the direction of the scattered photon and energy of 
the scattered electron. If the direction of the scattered elec- 
trons can be measured, the “event circle” becomes an “event 
arc”. These events are referred to as “tracked” events, and the 
length of the arc reflects the accuracy of the directional mea- 
surement of the scattered electron. However, measuring such 
events is difficult for VLAST because the scattered electrons 
are easily absorbed by the CsI low-energy detector. 

The STED has a total thickness equivalent to one radia- 
tion length, which allows for a 65% conversion rate of high- 
energy photons at normal incidence into electron-positron 
pairs. Fig. 2 illustrates the conversion positions. Fig. 3 pro- 
vides an orthogonal cut view of a 50 MeV gamma-ray event 
detected in the tracker, demonstrating the generation of an 
electron-positron pair in the CsI crystal and the gradual de- 
viation of the electron and positron from their original paths 
owing to multiple scattering. Recoil particles cause signal 
scattering. This inverted “V” signature is helpful in rejecting 
the much larger background of charged cosmic rays. 

The main design challenge of STED is to balance the elec- 
tron—positron pair conversion efficiency with the angular res- 
olution in the low-energy range. The angular resolution of 
the low-energy photons is significantly affected by multiple 
scattering, which depends on 1/E. Therefore, it is necessary 
to design thinner CsI conversion materials to reduce multiple 
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Fig. 2. The position distribution of 50 GeV photons converting into 
positron-electron pairs. From left to right, the peaks correspond to 
the ACD, STED, and HEIC, respectively. The cumulative distribu- 
tion function (CDF) of conversion probabilities is displayed by the 
y-axis on the right side of the figure. 
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Fig. 3. Illustration of the hits in the STED for a 50 MeV gamma-ray 
photon (x — z plane). The photon converts into an electron-positron 
pair in one of the CsI plates. 


scattering effects. However, this approach reduces the conver- 
sion efficiency of low-energy photons. To resolve this trade- 
off, CsI is divided into two types: thin layers are placed in 
the front section of the STED to ensure angular resolution for 
low-energy photons, and a thick layer is placed in the back 
to maintain a high conversion efficiency for high-energy pho- 
tons. For low-energy photons with energies of approximately 
100 MeV, the direction is mainly determined by the two lead- 
ing measurement points. The first two measurement points 
after the photon transition point are crucial in accurately re- 
constructing the direction of incident photons. To reduce the 
multiple scattering effects caused by the CsI plate and support 
material, the efficiency of each measurement layer should 
be close to 100%. Therefore, the spaced placement of CsI 
planes with silicon strip planes allows the precise measure- 
ment of positive and negative electron traces. The thick CsI 
plane placed at the bottom increases the efficiency of photon 
conversion at lower energies. This increases the statistics of 
low-energy photons. Studying the time-domain variation of 


photons is essential, despite the reduced angular resolution of 
photons. 


Another challenge is balancing the field of view (FoV) and 
angular resolution of the detector. For very-high-energy pho- 
tons, the impact of multiple scattering on the angular reso- 
lution becomes less critical. The ratio between the width of 
the silicon strip and the thickness of the silicon tracker detec- 
tor place the primarily limitation on the angular resolution. 
Ideally, finer silicon strips would yield improved angular res- 
olution. However, this results in an increased complexity of 
the fabrication process and higher electronic power consump- 
tion. By contrast, increasing the thickness of the STED leads 
to a reduction in the FoV and shifts the center point of the 
detector upward. STED was designed with these factors in 
mind. 


C. High Energy Imaging Calorimeter 


The HEIC has two main purposes. First, it measures the 
deposited energy of particles resulting from the interaction of 
incident photons. Second, to image the shower development 
profile which is used to reject the cosmic ray background and 
provide an estimate of the energy leakage fluctuations in the 
shower. To achieve high energy resolution, we adopted the 
design of the BGO calorimeter onboard the DAMPE [25]. 
The HEIC comprises four identical modular towers arranged 
in a 2x2 array; each tower has 416 BGO crystals divided 
into four layers, with each adjacent layer placed orthogonally. 
Each BGO crystal is 2.5 cm x 2.5 cm x 1.4 m in size. The 
total vertical depth of HEIC is 18 radiation lengths or 1 nu- 
clear action length, with oblique incidence events experienc- 
ing higher radiation lengths. The total effective detection area 
of HEIC is at last 2.8 m x 2.8 m. Photomultipliers at both 
ends read out the signals from the BGO bars. The different 
dynodes of the photomultiplier tubes are responsible for dif- 
ferent energy bands, enabling VLAST to encompass a wide 
energy range from MeV to TeV. The left/right light asymme- 
try provides a measure of the position of the energy deposit 
along the bar. Thus, each fired BGO crystal provides the x, y, 
and z coordinates of the shower, and the shower shape can be 
inscribed in three dimensions. The number of inscribed pix- 
els is directly related to the crystal dimensions. The shower 
axis also provides a rough track that can be used as a seed for 
the STED track reconstruction. Energy leakage is inevitable 
for grazing incidence and very high-energy photons, because 
the calorimeter has a limited radiation length. Leakage en- 
ergy can be estimated using a shower profile from beam tests 
or simulations. Taking the energy leakage correction into ac- 
count significantly enhances the energy resolution of these 
events. 


Ill. SIMULATION AND EVENT RECONSTRUCTION 
A. Simulation 


To validate and optimize the VLAST design concepts, we 
developed a simulation framework that includes incident par- 
ticle definition, detector geometry definition, a physical in- 
teraction list, digitization, and an event reconstruction analy- 
sis system. The framework is based on GEANT4, a publicly 
available toolkit developed in C++ used to simulate the inter- 
actions between particles and matter. It includes numerous 
physical models of particle-matter interactions and is widely 
used in various fields, such as high-energy physics, acceler- 
ator physics, space science, and medicine [49, 51-54]. The 
gamma-ray source is set up as a sphere with an energy range 
of 0.1 MeV to 1 TeV. 

The GEANT4 FTFP_BERT physics list was used to sim- 
ulate incident particles with energies above 10 MeV. For 0.1 
MeV-10 MeV gamma rays, we substituted electromagnetic 
physics with G4EmStandardPhysics_option4. This option in- 
troduces a more accurate model of Compton scattering and 
low-energy electromagnetic interactions, albeit at the cost of 
significantly increasing the computation time. We used the 
geometry module of GEANT4 to build the entire detector sys- 
tem, including the ACD, STED, and HEIC. The constructed 
VLAST geometry with a normal-incident 50 GeV gamma ray 
is shown in Fig. 4. The simulated output data of different sub- 
detectors were digitized with the typical electronic noise of 
their corresponding readout systems[55]. 


Fig. 4. Simulated shower particle tracks for a normal-incident 
50 GeV photon. The backsplash particles would produce self-veto 
in the ACD. 


B. Trigger design 


The preliminary trigger logic for scientific data collec- 
tion consists of five trigger engines: MeV-Gamma, GeV- 
Gamma, low-energy (LE), high-energy (HE), and Calibra- 
tion, as shown in Table 1. Each subdetector provides one or 
more trigger requests as detailed in the following list: 


e ACD: energy > 0.8 MeV, corresponding to 0.4 MIP 
(minimum ionization particle); 


CslI-hit | HEIC-hit | HEIC-HE | HEIC-LE | HEIC-MIP| STED | ACD 
MeV-Gamma 1 0 x x x x 0 
GeV-Gamma x 1 0 x x 1 0 

LE x 1 0 1 x 1 1 
HE x 1 1 x x x x 
Calibration x 1 0 x 1 1 1 


TABLE 1. Trigger logics of VLAST, including MeV-Gamma, GeV-Gamma, LE, HE and Calibration, and their corresponding requests (1: 


required; 0: excluded; x: either). 


e STED: three consecutive layers are on fire; 
e Cs]-hit: two layers are on fire; 
e HEIC-hit: > 5 MeV for the first or second layer; 


e HEIC-HE: > 5 MeV, > 500 MeV, > 500 MeV, 
> 500 MeV for the first 4 layers; 


e HEIC-LE: > 5 MeV, > 5 MeV, > 50 MeV, > 50 MeV 
for the first 4 layers; 


e HEIC-MIP: > 5 MeV for the Ist, 2nd, 7th, and 8th 
layers. 


The MeV-Gamma, GeV-Gamma, and HE trigger logics were 
combined to reserve gamma-ray events from sub-MeV to 
multi-TeV. The comprehensive trigger efficiencies were 45%, 
90%, and >95% at 1 MeV, 1 GeV and above 10 GeV, re- 
spectively. LE trigger logic (with a large prescale factor) was 
designed to reserve GeV cosmic rays for performance vali- 
dation in orbit. Calibration trigger logic was used for unit 
calibration of HEIC. 


C. Event reconstruction 


We utilized different algorithms to reconstruct the trajec- 
tories and energies of the pair production and Compton scat- 
tering events. The Kalman filter was used to reconstruct the 
trajectory of the pair production event after charge sharing. 
The Compton scattering formula was adopted to reconstruct 
the direction of photon incidence. The longitudinal develop- 
ment of an electromagnetic shower was used to reconstruct 
the energy. 

Before trajectory reconstruction, the process of charge 
sharing between the silicon strips should be considered. This 
is mainly due to diffusion during charge collection and ca- 
pacitive coupling. The drifting of electron-hole pairs along 
the magnetic field lines causes the size of the charge cloud 
to increase, whereas capacitive coupling occurs on the cou- 
pling strip between the two readout strips. The charge-sharing 
parameters originating from DAMPE were obtained using a 
beam test [25, 56-60]. The charge-sharing algorithm is as 
follows: 


E; = E; + 0.023- (E; + Bigs) + 0.0021 - (E;—4 + Ba); 
(2) 


if the strip is a readout strip and 


Ey =E; + 0.305 - (Ei—ı + Ei+1) + 0.062 - (Ei-3 + Ei+3) 
+ 0.012 - (Ei—5 + Ei+5) + 0.0024 - (Ei—-7 + Ei+7). 

(3) 
if the strip is a float strip. After the process of charge sharing, 
we combined the signals from neighboring silicon strips with 
a signal/noise ratio larger than four (19 keV) into a cluster, 
and then took the energy-weighted center as the position of 
the cluster. These clusters were used in the subsequent trajec- 
tory reconstruction analysis. 

In this study, the Kalman filtering algorithm was employed 
to reconstruct the trajectory of the pair production event, a 
technique widely used in particle physics experiments [61-— 
63]. When a gamma ray is converted into a pair in the CsI 
plane, the direction of the resulting electron/positron will ex- 
perience a shift owing to multiple scattering effects during 
propagation. Hence, understanding the impact of multiple 
scatterings on the trajectory-reconstruction process is crucial. 
The Kalman filtering algorithm assesses and compensates for 
multiple scattering and measurement errors. This algorithm 
consists of three primary processes: prediction, filtering, and 
smoothing. The track direction in the kth layer was used to 
predict the hit position on the (k + 1)th layer. Subsequently, 
the predicted hit position of the (k + 1)th layer was adjusted 
using the measured hit. The evolution of the state vector is 
given by: 

Xk = Fk-1Xk-1 + We-1, (4) 
where X;, is the state vector that incorporates the position and 
momentum information in the kth layer, Fķ—1ı is the propaga- 
tion of the trail from the (k — 1)th layer to the kth layer of the 
trail detector, and W;,_; is the random noise of the system. In 
dense media, tracking particles are subjected to random noise 
from multiple scattering, energy loss, and other physical pro- 
cesses that alter their trajectories. The trajectory offset caused 
by multiple scattering can be expressed as: 


13.6 
g = 


Bep 


where L is the thickness of CsI tile, L, is the radiation length 
of CsI material, Gc and p are the velocity and momentum 
(MeV) of the electron-positron pair, respectively. The mea- 
surement state vector is given by: 


L/L, [1 + 0.038 In(L/L,)], (5) 


Me = HkXpk + Ve, (6) 


where mz, is the quantity measured by the kth layer, H;, is the 
measurement matrix, and V; is the measurement error. Once 
a track had undergone a filtering process, it was smoothed. 
The trajectory parameters were further refined from bottom 
to top, in contrast to the filtering process. For further details, 
please refer to. [62, 64, 65]. 

The primary process for reconstructing the trajectory of 
MeV gamma-rays involves identifying Compton scattering 
events and determining the sequence of scattering points. The 
CsI hit found from top to bottom with an energy deposition 
greater than 100 keV and no readout signals (less than 19 
keV) from the silicon strips in the adjacent upper layer was 
considered the Compton scattering point. 

Using the same method, additional isolated scattering 
points can be identified, and their scattering sequence and the 
probability of photoelectric effects can be analyzed based on 
the magnitude of energy deposition. We traced the photon tra- 
jectory after Compton scattering by connecting the first and 
second interaction points with the highest probability. Ac- 
cording to the Compton scattering formula (1), reconstruc- 
tions of the angle of Compton scattering does not determine 
the specific direction of the primary photon, but only fixes a 
ring in the sky. Multiple rings from the same source can be 
overlaid to locate the source. 

For the electromagnetic shower induced by a high-energy 
photon, the overall deposited energy in CsI and BGO occu- 
pies most of the primary energy, despite a small fraction of 
energy loss. However, with an increase in incident energy, 
the longitudinal energy leakage becomes negligible, making 
the energy correction necessary to estimate the primary en- 
ergy. The longitudinal segmentation of CsI and BGO allows 
a fit of the longitudinal shower profile, which provides a good 
way to correct the longitudinal energy leakage [66]. The lon- 
gitudinal shower profile can be described using a gamma dis- 
tribution formula, expressed as: 

a-1. 8. Bt 
Og, OE Ee a 
t Tr(a) 


where t = x/ Xo represents the radiation length, T (œ) is the 
gamma function, œ and ĝ are shape and scale factors, respec- 
tively. The depth of the shower maximum depends on & and 
b as tmax = (a— 1)/2. The value of tmax is closely correlated 
to the energy leakage ratio, offering an effective approach for 
energy correction [66]. 


IV. EXPECTED PERFORMANCE OF THE VLAST 
A. Effective area 


The effective area reflects the detection efficiency of the in- 
strument. It is equivalent to the area of an ideal absorber that 
detects the same number of events as the real detector, con- 
sidering event selection and reconstruction, within the same 
time. The effective area is calculated as [67] 


Aer (E, 0, o) = Ageo(9, P)Edet (E, 0, P)Esa lE, 0, Q), (8) 


where Ageo(0, p) represents the geometric cross-sectional 
area of the instrument in a specific direction characterized by 
angles 0 and ¢, £qet( E, 0, ¢) is the detection efficiency, and 
Ese. (E, 0, Q) is the selection efficiency. The detailed selection 
conditions for pair events include the following: 1) the trigger 
condition for particular event types is satisfied (MeV-Gamma, 
GeV-Gamma, or HEIC-HE); 2) the trajectory in the STED 
can be reconstructed correctly; and 3) the deposited energy 
in the ACD block in the direction of the reconstructed trajec- 
tory is below 0.5 MeV. The selection conditions for Compton 
events include the following: 1) no ACD readout exceeding 
100 keV; 2) there are at least two CsI hits where the energy 
deposition at both points is greater than 100 keV and there 
is no readout signal (less than 19 keV) in the adjacent upper 
and lower layers of the silicon tracker detector; and 3) the en- 
ergy deposition in the bottom CsI layer is less than 50 keV (to 
ensure that the majority of photon energy is deposited in the 
STED). 


The effective area of VLAST is shown in Fig. 5. For 
Compton scattering events, the effective area is several thou- 
sand cm?, surpassing the effective area of the COMPTEL de- 
tectors (10-50 cm?) by one or two orders of magnitude. It 
also has a larger effective area than planned MeV detectors 
such as the ee ASTROGAM, AMEGO, GECCO, and XGIS- 
THESEUS detectors. Furthermore, the effective area for pair 
production events exceeds 4 m? above 1 GeV under normal 
incidence, which is four times larger than that of the Fermi- 
LAT detector. However, the effective area slightly decreases 
for 30° and 45° incidence angles for both Compton scatter- 
ing and pair production events. Fig. 6 shows the effective area 
as a function of the particle incidence angle at different ener- 
gies (50 GeV and 100 GeV). The effective area diminishes 
gradually with an increase in the angle of incidence, and the 
maximum off-axis incidence angle can reach 70°. The FoV 
of VLAST is 2.5 sr. 


The acceptance of a gamma-ray detector refers to the por- 
tion of incoming gamma rays that the detector is capable of 
capturing or detecting. A larger acceptance indicates that 
the detector can detect gamma rays from a wider portion of 
the sky, thereby increasing its sensitivity for gamma-ray de- 
tection. The primary design concept of VLAST is to sub- 
stantially enhance the acceptance of the detector to capture 
more photons and achieve more precise measurements with 
increased statistical significance. Therefore, the acceptance 
of VLAST must be sufficiently high. Acceptance is de- 
fined as the integral of the effective area over the solid angle 
G(E) = fa Act (E,9,¢)dQ. The acceptance of VLAST is 
shown in Fig. 5. For Compton scattering events, the max- 
imum value of acceptance is several m7sr. As the energy 
increases from MeV to GeV, the acceptance gradually rises 
from 1.6 m’sr to 12 m?sr. In the design of VLAST, CsI 
crystal plates replace the tungsten foils used in many other 
gamma-ray detectors to enhance pair conversion. This sub- 
stitution significantly improves the acceptance in the MeV 
energy range compared to Fermi-LAT. Given its substantial 
size, the acceptance of VLAST reaches 12 mĉsr in the en- 
ergy range above GeV. 
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Fig. 5. Left: effective areas of VLAST for photon events with normal incidence, 30°, and 45° incident angles, compared with those of 
AMEGO [12, 68] and Fermi-LAT [69]. Right: acceptance of the VLAST for gamma-rays of different energies, compared with results of 


GAMMA-400 [7] and Fermi-LAT [69]. 
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Fig. 6. Effective area for gamma rays with different incident angles. 


B.  Self-veto 


The detailed settings of the ACD need to be optimized, in- 
cluding the detection threshold and the size of the ACD. The 
detection threshold of the ACD must strike a balance between 
the detection efficiency of charged particles and suppression 
of backsplash effects. To improve the efficiency of detect- 
ing charged particles that exhibit a Landau distribution in the 
energy deposition within the ACD, the detection threshold 
should be lowered. However, to reduce the weak signals from 
recoil particles, the detection threshold should be set higher. 
The conflict between these two requirements is illustrated in 

Fig. 7. We used the false veto ratio of gamma rays to deter- 
mine the threshold. Gamma ray photons may be misidentified 
as charged particles and subsequently rejected because of the 
backsplash effect. Even if photons do not leave a signal when 
passing through the ACD, low-energy recoil photons can trig- 
ger the ACD. Consequently, the tile in the direction of the 
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Fig. 7. Energy deposition distribution of MIPs and backsplash effect 
in the ACD. The optimized threshold that balances the high detec- 
tion efficiency of charged primary particles and low contamination 
of backsplash is marked out in the plot. 


particle could exhibit a signal, leading to the event being in- 
correctly identified as a charged particle. The false veto ratio 
is defined as f = Nwrong /Ntotal, where Nwrong are the pho- 
tons misidentified as charged particles and Nota) is the total 
number of incident photons. The false veto ratio as a function 
of the energy is shown in the left panel of Fig. 8. The higher 
the photon energy, the greater the probability that it will be 
misclassified as a charged particle. The middle panel of Fig. 8 
presents the gamma-ray false veto ratios at different detection 
thresholds. For an ACD of size 20 cm x 20 cm, we deter- 
mined that the optimized threshold value was 0.5 MeV where 
the gamma false veto rate was less than 15%. The size of the 
ACD should also be set reasonably. Under ideal conditions, a 
smaller ACD block results in a better suppression of the back- 
splash effect. However, this also requires more readouts and 
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Fig. 8. False veto ratios as functions of incident energy (left), the ACD threshold (middle), and the ACD size (right). 


higher power consumption. Thus, to achieve the desired per- 
formance under limited conditions, tests were conducted to 
evaluate the suppression of the backsplash effect for different 
ACD sizes. The results of these tests are shown in the right 
panels of Figs. 8. An ACD tile size of 20 cm x 20 cm is rea- 
sonable when the false gamma veto rate is below 15%. Op- 
timization of the ACD thickness is also necessary. The ACD 
must be sufficiently thick to generate adequate fluorescence, 
ensuring the effective detection of charged particles. How- 
ever, it should not be excessively thick because this could lead 
to a significant number of photons interacting within the ACD 
being misidentified as charged particles. In this context, we 
set the thickness of the ACD to 1 cm, which is the same as that 
used for the Fermi-LAT. At this thickness, no more than 3% 
of the photons interact within the ACD. The interactions be- 
tween photons within the ACD are illustrated in Fig. 2. More- 
over, itis advisable to maintain wrapping materials around the 
ACD, such as micrometeoroid shields and thin thermal blan- 
kets. This precaution is important because charged cosmic 
rays can generate secondary photons that contribute to the lo- 
cal photon background upon interaction with such materials. 
This background can interfere with the accurate detection of 
gamma rays, thereby reducing the detection efficiency. 


C. Angular resolution 


Angular resolution is another critical parameter of VLAST. 
A better angular resolution not only provides more accurate 
source positioning, but also results in a sharper profile of the 
point source. This, in turn, leads to a smaller contribution 
from the background, thereby increasing the sensitivity for 
detecting fainter sources. The angular resolution is defined 
differently for different interaction processes, namely, Comp- 
ton scattering and pair production processes. For Compton 
events, the angular resolution is defined as the sigma value 
obtained from the Gaussian fit to the distribution of the min- 
imum angular distance between the nominal source position 
and the reconstructed event circle. For pair production events 
it is defined as the radius of the circle that includes 68% of the 
point-spread function (PSF). The PSF is determined using the 
angular distance between the source and the reconstructed di- 
rection. 


The angular resolution of VLAST is shown in Fig. 9. To 
obtain these results, the same selection conditions as those 
described in Sec. IV A for calculating the effective area are 
used. For Compton events, the angular resolution is within 
the range of 4° to 8°, which is comparable to AMEGO. 
Here, we present only the angular resolution of the untracked 
events. For pair production events, the angular resolution im- 
proves with increasing photon energy owing to the reduced 
impact of multiple scattering effects at higher energies. No- 
tably, at 10 GeV, the angular resolution reaches 0.2 degrees. 
However, in the energy range of 10 to 10° MeV, the angular 
resolution is slightly worse than that of Fermi-LAT. In con- 
trast, in the energy range above GeV, the angular resolution 
surpasses that of Fermi-LAT. This improvement can be at- 
tributed to our use of analog readout, whereas the Fermi-LAT 
employs a digital readout with silicon strips of approximately 
the same width. 
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Fig. 9. Angular resolution (68% containment) of VLAST for photon 
events with normal incidence, 30° and 45° off-axis angles, com- 
pared with results of AMEGO [12, 68] and Fermi-LAT [69]. 


As mentioned earlier, we implemented a spaced readout 
scheme for the silicon strip readout to reduce power con- 
sumption. However, VLAST requires a large number of elec- 
tronic readouts, which results in substantial power consump- 
tion. To address this issue, we explored the use of digital 
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Fig. 10. Angular resolution of VLAST for normal incident photons. The left panel shows the comparison between digital and analog readouts 
for a 120 um strip width, and the right panel compares different widths of silicon strips for analog readout. 
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Fig. 11. Angular resolution (left) and effective area (right) of VLAST for different configurations of CsI layers for Compton events. Normal 


incidence and analog readout are adopted in the simulation. 


readout, which records only the hit strips that exceed a sig- 
nal threshold of 30 above the baseline noise as “1”, while 
unhit strips are recorded as “0”. This implementation simpli- 
fies the electronics significantly and reduces the power con- 
sumption. We defined a cluster as the region of the hit silicon 
strips. For the digital readout, during cluster position recon- 
struction, the cluster position was set to the geometric cen- 
ter, whereas for analog readout, the energy-weighted average 
position was used for cluster reconstruction. Subsequently, 
track reconstruction was performed using Kalman filtering. 
A comparison of analog and digital readouts is shown in the 
left panel of Fig. 10. In the lower-energy range, the domi- 
nant error arises from multiple scattering, resulting in com- 
parable angular resolutions for both analog and digital read- 
outs. However, in the higher-energy range, measurement er- 
rors take precedence, and the geometric center of the cluster 
exhibits a significantly larger error than the energy-weighted 
center. This difference increases with energy, indicating that 
analog readouts provide better angular resolution than digital 


readouts. 


We investigated the effect of different silicon strip widths 
on the PSF of VLAST for pair production events. Thinner 
strips offer a higher position measurement accuracy and bet- 
ter angular resolution. However, they require more electronic 
readouts, resulting in higher power consumption. Therefore, 
a balance between the angular resolution and power con- 
sumption is required. The tests were conducted using three 
widths: 80 um, 120 um, and 160 um. The results are shown 
in the right-hand panel of Fig. 10. Notably, there was no 
significant difference in the angular resolution in the lower- 
energy band among the different widths. This is because the 
uncertainty in the angular resolution is primarily attributed to 
multiple scattering effects within this energy range. However, 
as multiple scattering effects decrease at higher energy levels, 
the measurement uncertainty becomes more prominent, and 
finer silicon strips lead to better resolution. As a benchmark, 
we chose a silicon strip with a width of 120 um, which pro- 
vided an angular resolution of 0.05° at 50 GeV. 


In this study, we investigated the effect of different CsI 
layers on the angular resolution. To simplify the mechani- 
cal structure and ease the manufacturing and assembly pro- 
cesses, the STED favors thicker CsI plates with fewer layers. 
However, thicker CsI plates degrade the angular resolution 
owing to the larger deflection angle resulting in multiple scat- 
tering effects. Therefore, the thickness of the CsI plate must 
be set appropriately. We tested three different configurations 
with 4, 8, or 16 layers of CsI plates interleaved in the STED 
while maintaining the same total radiation length to ensure 
the same efficiency of pair production. The effects of differ- 
ent layer settings on Compton-scattering events are shown in 
the left panel of Fig. 11. Different layer settings had oppo- 
site effects on the angular resolution and effective area. The 
angular resolution decreased with increasing number of CsI 
layers. Because the four layers had the thickest CsI, the mea- 
sured energy was the most accurate, and the corresponding 
reconstruction direction was the most accurate according to 
the Compton scattering equation. However, the effective area 
increased with an increase in the number of CsI layers. The 4- 
layer CsI configuration exhibited the smallest effective area. 
For pair production events, the effects of different CsI layer 
configurations on the angular resolution are shown in Fig. 12. 
The angular resolution was not significantly different in the 
low-energy band. The 16-layer configuration offers the best 
angular resolution at several GeV because of reduced multi- 
ple scattering in the thinner CsI. However, at several hundred 
GeV, the difference becomes negligible as multiple scatter- 
ing effects diminish, whereas the measurement errors become 
more significant when the same pitch of the silicon strip is 
used. Different layer configurations had the same effective 
area because the total radiation length was the same. Based 
on comprehensive considerations, we determined that setting 
the number of CsI layers to 8 yielded the most suitable bal- 
ance. 
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Fig. 12. Angular resolution of VLAST for different configurations 
of CsI layers for pair events. Normal incidence and analog readout 
are adopted in the simulation. 
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D. Energy resolution 


Energy resolution is another critical aspect of the VLAST. 
One of the most interesting scientific objectives of the 
VLAST is the line structure in the spectra produced from 
astrophysical phenomena, such as de-excitation nuclear 
gamma-ray line emission from low-energy cosmic rays and 
gamma-ray lines from the annihilation of dark matter parti- 
cles [70-72]. The ability to detect the line structure in the 
gamma ray spectrum relies significantly on the energy reso- 
lution. Here, we analyzed the energy resolution of the STED 
and HEIC using the above-mentioned configuration. Firstly, 
for pair events, the trigger conditions are MeV gamma, GeV 
gamma, and HEIC-HE. Second, in STED, Kalman filtering 
can reconstruct the photon trajectory. Third, most electro- 
magnetic showers are contained in the HIEC. The energy de- 
position was obtained from the simulation, considering the 
statistical fluctuation error and electronic noise (based on 
DAMPE [73-76]) was employed for energy resolution anal- 
ysis. The energy resolution was obtained after energy correc- 
tion for high-energy events with more energy leakage. For 
Compton events, the energy was measured exclusively by the 
CsI calorimeter, and the uncertainty of the measurements was 
primarily considered in four aspects: quantum fluctuations 
of CsI scintillation light yield, transmission efficiency of the 
WLS fibers, photon detection efficiency and electronic noise 
of SiPMs, and the accuracy of energy calibration. The energy 
measurement process was simulated by randomly sampling 
each energy deposition of CsI using specific parameters ini- 
tially estimated from previous experiments. The energy res- 
olution was obtained using a Gaussian fit of the proportion 
of the sample results to the simulated real energy. The energy 
resolution is shown in Fig. 13. For the pair production event, a 
peak at approximately 100 MeV is present because the depo- 
sition body shifts from a low-energy calorimeter to the HEIC. 
As the number of energy particles entering the calorimeter in- 
creases, the improving energy resolution reaches 2% at a few 
tens of GeV. Subsequently, the resolution becomes weak ow- 
ing to energy leakage. The energy resolution of the 30° and 
45° incidence events were better than those of the normal- 
incidence events because the radiation length of the path was 
longer. Compared with Fermi-LAT, which has an approxi- 
mate radiation length of 8.6, VLAST exhibits better energy 
resolution across the entire energy band owing to its deeper 
calorimeter, except at approximately 100 MeV. 


E. Electron and proton discrimination 


For effective photon detection, VLAST requires a strong 
capability to reject the charged-particle background. As pre- 
viously mentioned, the proton flux is five orders of magnitude 
higher than that of the photons, making it difficult to detect 
gamma rays in the presence of such a large background of 
charged particles. In addition to the ACD, the HEIC plays a 
crucial role in background rejection by distinguishing protons 
from photons based on the difference between hadronic and 
electromagnetic showers. To quantify the e-p discrimination 
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Fig. 13. Energy resolution of VLAST for Compton and pair events 
with normal incidence, 30° and 45° incident angles, compare with 
results of AMEGO [12, 68] and Fermi-LAT [69]. 


ability of the VLAST, we used the ¢ parameter, as in DAMPE 
[77], which is defined as 


¢ = Fast X (©iRMS;/mm)*/(8 x 10°), (9) 


where Fj,5; is the ratio of energy deposition in the last layer to 
the total energy deposition in the HEIC. The energy deposi- 
tion in the last layer reflects the differences between hadronic 
and electromagnetic showers most significantly because BGO 
crystals have a large nuclear interaction length-to-radiation 
length ratio. RMS; is the root-mean-square value of the en- 
ergy deposited hit position in the zth layer, expressed as: 


0 (25,4 — Lei)? Eji 
ROS. = yj Mara oa) Bie (10) 


455 


where x; ; and E; ; are the positions of the hits and the de- 
posited energy of the j-th bar in the i-th layer, respectively, 
and x, ; is the center coordinate of the shower in the 7-th layer. 
This reflects the lateral development of the shower, which 
is primarily caused by the propagation of secondary parti- 
cles. The secondary particles produced by hadronic showers 
can propagate farther than those produced by electromagnetic 
showers [78], which results in larger values of RMS;. 

The left panel of Fig. 14 shows the scattering distribution 
of RMS = >>, RMS, and Fj,s¢ for protons (orange) and elec- 
trons (blue), respectively. Both the electron and proton events 
have energy depositions in the HEIC range of 30-100 GeV. 
Electrons and protons were clearly divided into two parts. 
The distributions of ¢ parameters for the electron and proton 
samples are shown in the right panel of Fig. 14. A ¢ value 
of 18 was determined to suppress proton contamination to 
<0.1% while maintaining 90% electron detection efficiency. 
By choosing proper cuts in the ¢ parameter, we can obtain 
the rejection fraction of the protons while maintaining 90% 
of the electrons, as shown in Fig. 15. At relatively low ener- 
gies (several GeV), the proton rejection fraction was approx- 
imately 0.99 owing to the small size of the electromagnetic 
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cascades. With an increase in energy, the rejection fraction in- 
creases effectively, reaching 0.999 at ~ 20 GeV. Optimization 
of the e-p discrimination algorithm, for example, by means of 
machine learning, can further improve the background rejec- 
tion capability [79]. Combined with the background rejection 
fraction of ~ 0.999 from the ACD (design requirement) and 
~ 0.9 from the STED, and based on the experience of Fermi- 
LAT [5], the total background rejection fraction of VLAST 
can reach (1 — 1076) — (1 — 1077). Thus, VLAST has an 
excellent ability to detect gamma rays, even if the cosmic-ray 
background flux is 10° higher than that of the photons. 

However, VLAST was not optimized for electron detec- 
tion. The vertical thickness was approximately 18 radiation 
lengths, resulting in a relatively large leakage of electron 
events with energies above | TeV. However, given the large 
area of the detector, events with large incident angles can be 
selected to effectively increase the slant thickness of the de- 
tector. As a reference, for incident angles > 55°, the thick- 
ness is approximately 32 radiation lengths, which is compa- 
rable to that of DAMPE. If we select events with incident 
angles between 55° and 75°, the acceptance (for the HEIC- 
HE trigger) is approximately 3.5 m?sr. Because DAMPE can 
measure the electron spectrum up to approximately 15 TeV, 
VLAST can extend the detectable upper energy bound by at 
least a factor of two, even if the spectrum is as soft as E~4 
[77]. However, the dynamic range of electronics must also be 
expanded to detect such high-energy events. 


F. Detectablity of gamma-ray transients 


VLAST has a very good capability for exploring the y-ray 
sky with an unprecedented sensitivity over a wide range of en- 
ergies [21]. Considering the detection of burst-like transients 
as an example, we briefly discuss the potential of VLAST for 
such phenomena. Assuming a physical sample of gamma ray 
bursts (GRBs) with parameterized distributions of redshifts, 
luminosity, and spectra, the expected detectable quantities of 
GRBs can be calculated by a given detector using parame- 
ters such as the effective area, energy band coverage, PSF, 
and FoV [81-83]. The intrinsic parameters were obtained by 
comparing the simulated results with the observed sample us- 
ing various detectors. The detectability can then be calculated 
for new detectors such as VLAST. Fig. 16 shows the expected 
yearly rate for different fluences using VLAST and compares 
the expected rate by Fermi-LAT [81] and the Fermi-LAT de- 
tection rate [80] above 100 MeV. The total annual rate is ap- 
proximately 43.4 for VLAST, which is approximately three 
times higher than that of Fermi-LAT (14.5 per year). 


V. SUMMARY 


To validate and optimize the design of VLAST, we simu- 
lated and analyzed the performance parameters of VLAST 
and optimized some detailed settings such as the effective 
area, angular resolution, energy resolution, threshold and size 
of ACD, and width of the silicon strip. WLAST directly in- 
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Fig. 15. The proton rejection fraction as a function of the energy 
deposition in HEIC when retaining 90% efficiency of electrons. 


creases the gamma-ray detection capability by increasing the 
effective area while replacing the conventional tungsten plate 
used for electron pair conversion with CsI to detect Comp- 
ton scattering events in the MeV band. This design allows 
the VLAST to have an effective area of 4 m?, which is larger 
than the previous ones and smaller than the APT, but with 
better energy and angular resolution. The outstanding perfor- 
mance of VLAST can identify the uncertified point sources 
of Fermi-LAT, further distinguishing whether the gamma-ray 
excess at the center of the Milky Way originates from dark 
matter, whether the Fermi bubble is of leptonic or hadronic 
origin, and whether it fills in the gaps in the MeV band of 
the extragalactic background light. The design and valida- 
tion of VLAST are continue being a subject of research. In 
future work, we consider incorporating additional subdetec- 


tors such as time-of-flight and neutron detectors to improve 
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Fig. 16. Expected detectable rate of GRBs with different fluences by 
VLAST (red line), compared with the detected sample [80] and the 
expectation [81] by Fermi-LAT. 


the rejection capability for complex backgrounds. Time-of- 
flight detectors can reduce the backsplash effect by measur- 
ing the time of recoil photons, whereas neutron detectors can 
improve the background rejection capability of the VLAST 
because hadron interactions produce large numbers of neu- 
trons compared to electromagnetic interactions. More de- 
tailed event analysis algorithms are also under development to 
obtain a more accurate performance of VLAST. A prototype 
for validating the principles of the VLAST design is currently 
under development. VLAST is expected to play a crucial role 
in gamma-ray astronomy. 
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